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Summary.  This  research  program  involves  mathematical  and  statistical 
techniques  for  using  models  to  analyze  experimental  data,  to  aid  in  experi¬ 
mental  design,  and  to  optimize  system  designs.  The  applications  of  interest 
include  subsurface  contaminant  transport,  and  flow  measurement  and  opti¬ 
mal  design  in  aerospace  testing  experiments. 

Fate  and  Transport  Research.  Research  continues  on  efficient  modeling 
and  parameter  estimation  in  heterogeneous  media.  We  have  developed  new 
models  based  on  point  processes,  which  allow  for  “blocky,”  piecewise  constant 
heterogeneities,  and  we  have  tested  these  methods  against  tradition  kriging 
estimates  as  well  as  our  regularized  least  squares  approach.  We  have  also 
developed  transport  models  based  on  the  statistical  characterizations,  and  we 
have  developed  numerical  and  statistical  methods  for  analyzing  unsaturated 
flow  problems  for  heterogeneous  media.  Michael  Jeffris  completed  his  Ph.D. 
thesis  on  simulation  of  unsaturated  flow  in  heterogeneous  media. 

Inverse  Interferometry.  In  joint  work  with  Dr.  Steve  Keeling  of  AEDC, 
we  have  developed  and  tested  numerical  algorithms  for  estimating  flowfield 
density  from  interferogram  data.  The  algorithms  involve  a  minimum  residual 
cost  functional  penalized  with  a  total  variation  regularization  term.  We  have 
developed  techniques  for  distributing  the  regularization  level  in  an  adaptive 
manner,  to  capture  better  piecewise  smooth  density  profiles. 

Turbulence  Models  and  Optical  Propagation.  The  Airborne  Laser  (ABL) 
program  involves  the  propagation  of  laser  light  over  long  distances  in  the 
atmosphere,  and  during  its  travel  the  light  encounters  significant  turbulence. 
The  effects  of  the  turbulence,  which  manifests  itself  through  index  of  refrac¬ 
tion  variations,  are  to  spread  the  light  beam,  thereby  reducing  its  effective¬ 
ness.  A  major  challenge  for  this  program  is  the  statistical  characterization 
of  the  turbulent  atmosphere,  so  that  accurate  propagation  simulations  can 
be  made  to  assess  performance  of  system  designs.  In  collaboration  with  Dr. 
Don  Washburn  of  PL/LIA,  we  have  transitioned  our  statistical  methods  for 
site  characterization  in  contaminant  transport  applications  into  techniques 
for  understanding  the  turbulence  data  obtained  by  Phillips  Lab  scientists. 
These  methods  allow  us  to  extract  regions  of  stationarity  from  non-stationary 
time  series  data. 

Accomplishment  s . 
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Subsurface  contaminant  transport.  In  this  research  area,  we  have  made 
progress  in  several  modeling  and  parameter  estimation  problems  involving 
transport  in  saturated  heterogeneous  media. 

Most  of  our  effort  has  focused  on  a  type  of  random  process  known  as  a 
modulated  Boolean  point  process,  or  MBPP  for  short,  which  provides  a  sim¬ 
ple  statistical  model  of  “blocky”  functions  (i.e.,  functions  which  are  piecewise 
smooth,  with  fairly  complex  discontinuities).  These  processes  are  modeled 
statistically  using  two  probability  distributions:  one  for  the  number  of  blocks 
and  the  other  for  the  size  of  the  blocks. 

A  drawback  of  the  standard  second-order  stationary  modeling  used  for 
conductivities  is  that  the  realizations  do  not  typically  “look  like  the  true 
fields.”  One  expects  some  small  scale  variation  which  resembles  a  second- 
order  stationary  process,  with  some  larger  scale  variations  that  are  “blockier.” 
A  modeling  approach  that  provides  such  a  structure  is  the  modulated  point 
process. 

The  modeling  begins  with  a  (random)  number  of  points  in  the  region  of 
interest,  points  which  represent  centers  of  the  “blocks.”  The  next  step  is  the 
block  size  and  shape.  The  mathematical  formulation  of  this  idea  is  a  function 

i= 1  r* 

in  which  N  is  &  discrete  random  variable  denoting  the  number  of  blocks, 
Ki  denotes  the  “amplitude”  of  the  block,  &  denotes  the  random  block  cen¬ 
ters,  Ti  denotes  the  block  sizes,  and  the  (deterministic)  function  <f  gives  the 
block  shape.  For  example,  <f>  can  be  a  Gaussian  function,  yielding  smoothly 
shaped  blocks;  alternatively,  (f)  can  be  a  box-car  function,  producing  blocks 
with  sharp  edges.  The  function  (j>  can  also  contain  anisotropies  (imagine  a 
Gaussian  function  with  a  non-diagonal  covariance  in  it).  This  type  of  process 
is  (when  all  the  Ki  values  axe  one)  the  standard  Boolean  point  process.  The 
modulated  process  arises  from  putting  a  Boolean  point  process  “on  top  of’ 
the  original  blocks. 

The  mathematical  procedure  for  building  this  model  is  as  follows.  For 
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each  center  and  size  above,  we  have 


«(»)-* 

j=i  'u 

with  fa  given  by 

Then  the  modulated  process  is  given  by 

/(*)  =  Ey<W- 

i=l 

In  this  construction,  the  new  centers  £y  and  radii  ry  have  a  distribution 
dependent  on  the  values  &  and  ry  for  each  i.  Realizations  of  these  processes 
have  the  appearance  of  nonstationarity,  but  in  fact  the  process  is  a  stationary 
process.  Since  it  is  often  mentioned  in  the  literature  that  field  sites  do  not 
appear  to  be  stationary,  we  feel  that  this  modeling  tool  may  be  particularly 
appropriate  for  hydraulic  conductivity  modeling.  Our  results  for  the  NAS 
site  at  Columbus  AFB  are  quite  encouraging. 

We  should  also  remark  that,  while  the  realizations  of  this  type  of  process 
are  dramatically  different  from  standard  second-order  stationary  processes, 
the  covariances  can  be  quite  similar.  This  fact  has  implications  for  the  ‘self- 
similarity’  modeling  that  is  currently  being  proposed  to  explain  the  scale 
dependence  of  dispersion  parameters. 

The  estimation  problem  involves  estimating  the  values  of  the  random 
quantities  N,  the  block  radii  and  centers,  and  the  intensities.  We  have  de¬ 
veloped  a  Markov  Chain  Monte  Carlo  algorithm  to  generate  distributed  es¬ 
timates  of  the  hydraulic  conductivity  field  from  observations  of  the  field  at 
discrete  points.  The  method  works  as  follows.  One  initiates  a  random  walk 
modeling  the  center  location  of  the  modulator  function  (i.e.,  the  &),  with  a 
fixed  (and  large)  radius  (i.e.,  the  r*).  The  random  walk  is  like  a  simulated  an¬ 
nealing  search,  directed  toward  regions  with  datapoints  having  similar  values 
for  the  conductivity.  Once  this  point  is  established,  we  begin  another  ran¬ 
dom  walk,  and  we  keep  adding  modulator  sets  until  we  obtain  insufficient 
improvement  in  the  fit.  Then  we  reduce  the  radius  size  and  begin  again. 
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We  continue  this  procedure  until  the  radii  become  too  small.  Of  course,  “too 
small”  and  “insufficient  improvement”  mean  that  the  user  must  specify  these 
tolerances.  The  final  stage  of  this  effort  is  to  devise  robust  and  reliable  means 
of  assisting  the  user  in  specifying  these  toerances. 

One  interesting  point  is  that  it  is  a  straightforward  matter,  based  on 
the  particular  choices  of  probability  distributions  for  the  random  quantities 
in  the  model,  to  obtain  stationary  or  non-stationary  behavior,  as  well  as 
the  self-similar  or  ‘fractal’  behavior  proposed  by  Neumann  and  others  as  an 
explanantion  of  the  scale-effect  in  dispersion  modeling. 

In  our  contaminant  transport  modeling,  major  effort  has  continued  on 
modeling  hydrodynamic  dispersion.  A  serious  drawback  in  the  application 
of  standard  transport  models  to  field  scale  problems  is  the  determination  of 
the  dispersivity  parameters.  The  basic  model  for  a  conservative  tracer  in 
groundwater  is  given  by 


net  +  V  •  ( vc )  =  V  •  (D'Ve), 

in  which  c  =  c(t,  x)  is  the  contaminant  concentration,  n  is  the  porosity,  v  is 
the  groundwater  velocity,  and  D  is  the  dispersion  matrix,  given  by 

D  =  d0I  +  (Qz-  +  aT\v\I, 

M 

in  which  d0  is  the  molecular  diffusion  coefficient,  and  Ol  and  Or  are  the 
longitudinal  and  transverse  dispersivities.  The  difficulty  in  applying  this 
model  is  that  the  dispersivity  parameters  appear  to  depend  on  the  scale  of 
observation:  in  field  scale  problems,  estimated  dispersivity  increases  as  the 
size  of  the  problem  increases. 

The  approach  we  have  taken  to  investigate  this  problem  is  a  stochastic 
one:  modeling  the  variation  in  velocity  directly  and  determining  a  transport 
model  which  includes  a  collection  of  velocities.  It  is  the  variation  in  velocities 
at  multiple  scales  which  is  believed  to  produce  the  plume  spreading  which  is 
(hoped  to  be)  captured  by  dispersion  models. 

The  basic  idea  behind  the  above  advection-dispersion  equation  is  that  the 
variation  in  the  velocity  field  can  be  modeled  as  an  increase  in  the  (Fickian) 
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diffusion.  Theoretical  justifications  for  such  a  model  rely  on  large  time  and 
small  variance  asymptotics,  which  are  in  general  not  applicable  for  many 
heterogeneous  media.  Our  investigations  are  along  two  distinct  but  related 
lines.  One  is  related  to  the  point  process  methods  discussed  above;  the  other 
is  a  direct  Monte  Carlo  computation. 

The  point  process  model  begins  by  imagining  the  porous  medium  to  be 
a  sequence  of  well- mixed  reactors  (cells)  connected  in  series.  The  modeling 
assumption  that  we  make  is  that  the  time  spent  in  each  cell  is  random.  These 
random  times  are  assumed  independent  and  identically  distributed. 

We  denote  by  Ax  the  length  of  each  cell.  The  time  spent  in  the  ith  cell 
is  Tj.  We  denote  by  fxt  and  of  the  mean  and  variance  of  these  random  times. 
Since  most  physical  models  of  porous  media  treat  the  velocity  (directly,  or  in¬ 
directly,  through  the  hydraulic  conductivity)  as  the  “basic”  random  quantity, 
we  consider  random  velocities  Vi,  related  to  these  times  by  the  relationship 
Ax  =  VjTi .  We  denote  the  velocity  mean  and  variance  by  jiv  and  of,  re¬ 
spectively.  In  general,  one  must  know  the  probability  distribution  of  V*  to 
determine  the  mean  and  variance  of  T*,  since  the  relationship  between  the 
two  is  nonlinear.  Below  we  shall  give  an  example  of  how  to  relate  the  time 
and  velocity  parameters  for  lognormally  distributed  conductivity. 

The  time  required  to  pass  through  the  first  n  cells  is  given  by  Zn  = 
SiLi  Alternatively,  we  can  determine  the  cell  containing  the  particle  at 
time  t  by  M(t )  =  mi{n:Zn  >  t}.  The  distance  traveled  by  the  particle  in 
this  system  is  then  given  by  X(t)  =  AxM(t).  The  probability  density  of 
the  particle  position  is  the  quantity  related  to  the  solution  of  the  advection 
dispersion  equation. 


The  central  limit  theorem  for  renewal  processes  indicates  that 

M(t)  -  t/nt 


y/ott/tf 


where  the  convergence  is  convergence  in  distribution.  Of  course,  to  relate 
this  result  to  the  travel  distance,  X(t ),  we  simply  multiply  and  divide  by  Ax 
to  obtain 

Xj^-A 
\J  Ax2  o\tjn\ 
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The  asymptotic  mean  term  in  this  model  is  given  by 


E(X(t)) 


A  xt  _  Axt 


and  the  asymptotic  variance  is 

Vt 

We  point  out  here  that  the  asymptotic  mean  contains  the  term  l/E(l/Vi) 
which  is  of  course  different  from  E(Vi).  Our  numerical  studies  indicate  that 
this  parameterization  compares  much  better  to  Monte  Carlo  simulations  that 
the  commonly  used  approximations  based  on  particle  tracking  in  a  second 
order  stationary  velocity  field. 

Another  interesting  line  of  inquiry  has  been  in  the  study  of  probability 
models  for  the  velocity  that  have  heavy  tails.  In  fact,  using  results  from  the 
modeling  of  network  traffic,  the  anomalous  dispersion  described  by  Dagan 
and  Neumann  can  be  understood  in  that  context.  The  basic  model  of  the 
sequence  of  mixing  regions  remains  the  same,  but  instead  of  using  only  first 
and  second  statistical  moments,  we  assume  a  probability  distribution  whose 
tailhas  a  slow  decay  rate.  In  this  circumstance,  the  central  limit  theorem 
produces  a  different  limit,  whose  dispersion  coefficient  matches  field  data 
much  more  accurately. 

As  noted  above,  in  order  to  determine  these  parameters  from  the  phys¬ 
ical  parameters  (i.e.,  the  hydraulic  conductivity’s  statistical  properties),  we 
must  know  the  probability  distribution  of  the  velocity  field.  We  have  per¬ 
formed  these  computations  using  lognormally  distributed  conductivity  pro¬ 
cesses  (which  is  the  most  common  model  in  the  hydrogeology  literature),  and 
we  have  seen  very  good  agreement  with  the  Monte  Carlo  simulations  (which 
are  discussed  below). 

Our  second  approach  is  to  model  the  plume  as  being  a  superposition  of 
many  plumes,  each  following  a  different  velocity  field.  These  different  ve¬ 
locity  fields  are  modeled  as  random  fields  arising  ^from  random  field  models 
of  the  conductivity.  The  plume  is  then  determined  by  summing  up  these 
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“sub-plumes.”  The  computations  involve  the  Monte  Carlo  simulation  of 
the  hydraulic  conductivities,  which  through  Darcy’s  law  generate  velocity 
fields.  The  advection-diffusion  equation  (i.e.,  D  =  dol,  no  dispersion  terms) 
is  solved  using  these  random  velocity  fields  instead  of  an  average  velocity, 
and  the  “sub-plumes”  are  superposed  to  obtain  the  plume  behavior.  We 
have  proved  convergence  theorems  for  these  numerical  approximations  to 
capture  the  mean  behavior  of  the  plume,  and  we  have  conducted  extensive 
numerical  simulations  to  compare  with  simpler  approximate  models  found 
in  the  literature.  Much  of  this  work  can  be  found  in  our  recently  submitted 
paper  (in  the  publications  list  below).  An  advantage  to  this  approach  (which 
we  are  currently  investigating)  is  that  one  can  use  a  wide  variety  of  random 
process  models  for  the  conductivity.  We  are  in  the  process  of  comparing 
the  MBPP  and  second-order  stationary  models,  from  the  point  of  view  of 
transport. 

Inverse  Interferometry.  In  this  research  effort,  a  collaboration  with  Dr. 
Steve  Keeling  of  Sverdrup  Technology,  Inc.,  Computation  and  Modeling 
Branch,  AEDC,  Arnold  Air  Force  Base,  we  have  examined  some  theoretical 
issues  associated  with  the  use  of  total  variation  based  image  reconstruction. 
Our  investigations  are  motivated  by  problems  of  inverse  interferometry,  in 
which  laser  light  phase  shifts  are  used  to  reconstruct  medium  density  pro¬ 
files  in  flow  field  sensing.  The  reconstruction  problem  is  posed  as  a  residual 
minimization  with  total  variation  regularization  applied  to  handle  the  inher¬ 
ent  ill-posedness.  We  consider  numerical  approximations  of  these  penalized 
minimal  residual  problems,  and  analyze  some  approximation  strategies  and 
their  properties.  The  standard  definition  of  total  variation  leads  to  inconsis¬ 
tent  approximations,  with  piecewise  constant  basis  functions,  so  we  consider 
alternative  definitions,  which  preserve  the  needed  compactness  and  produce 
convergent  approximations. 

The  general  problem  to  be  considered  involves  a  theoretical  image  /  €  X, 
an  observed  image  g  e  Z,  and  an  observation  operator  $  :  X  -»  Z.  For 
instance,  X  =  Ll( F ),  where  F  is  a  bounded  flow-field  domain  in  Rd,  and 
Z  =  Ll(P),  where  P  is  a  bounded  measurement  domain  in  R2.  Alternatively, 
Z  may  be  a  finite  dimensional  space  of  pixel  values.  A  typical  problem 
involves  solving  for  /  in  $/  =  g  with  observed  data  g  €.  Z.  The  problem 
is  treated  by  minimizing  some  measure  of  the  residual,  $/  —  g,  subject  to 
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smoothing  constraints  which  de-noise  the  image,  /.  In  particular,  we  have 
studied  theoretically  and  computationally  penalized  cost  functionals  of  the 
form 

Mf)  =  W  -  g\\»  +  1&V1V), 

in  which  /  represents  the  difference  between  the  “far  field”  fluid  density  and 
the  density  in  the  region  of  interest,  g  represents  the  interferogram,  and  the 
forward  operator,  $,  is  the  interferometry  operator 

«(/)( [x,y)  =  Gf*-j&Xirdr, 

with  G  being  a  known  constant.  This  operator  is  derived  assuming  axisym- 
metric  flow. 

We  remark  here  that  the  use  of  the  Ll  residual  criterion  arises  for  sev¬ 
eral  reasons.  One  is  that  L 1  is  a  natural  space  for  flow- field  densities,  and 
the  operator  $  lies  in  C(Ll(F),Ll(P)).  Another  is  that  the  effectiveness  of 
smoothing  via  the  total  variation  requires  BV  to  be  compactly  imbedded  in 
X,  a  condition  met  by  X  =  L1  but  not  X  =  L2  in  dimension  d>  2.  Finally, 
the  Ll  criterion  is  typically  more  robust  in  a  statistical  sense  than  the  L2 
criterion. 

The  notation  TVi  indicates  that  our  definition  of  the  total  variation 
is  not  the  traditional  one.  The  details  of  the  definition  are  somewhat  te¬ 
dious;  it  is  most  easily  explained  for  smooth  functions,  for  which  TVi.  (/)  = 
/F  |V/(x)|i  dx  where  |  •  |i  denotes  the  1-norm  in  R2.  The  use  of  this  def¬ 
inition  of  total  variation  is  needed  in  order  to  devise  consistent  numerical 
approximations  for  piecewise  constant  functions.  In  order  to  make  the  cost 
functional  differentiable,  we  consider  the  approximation 

U(f)  =  Jf  V/72  +  |V/|2dx 

and  in  our  paper  we  prove  that  minimizers  of  the  cost  penalized  with  U 
converge  to  those  with  TVi  as  the  penalty,  as  7  — >  0. 

The  focus  during  this  reporting  period  has  been  on  advancing  the  numer¬ 
ical  work  and  improving  the  speed  and  accuracy  of  the  code.  One  particular 
area  of  effort  has  been  in  distributed  regularization.  Distributed  regulariza¬ 
tion  is  a  technique  which  uses  a  spatially  varying  regularization  parameter 
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to  weight  the  image  differently  in  regions  of  differing  smoothness.  We  are 
currently  investigating  the  use  of  a  penalty  term  of  the  form 


in  which  /0  is  a  “lagged”  estimate  of  the  density,  and  a;  is  a  weight.  The 
function  fi/yfl  +  w2|V/o|2  plays  the  role  of  the  regularization  parameter,  be¬ 
coming  nearly  zero  in  regions  of  large  gradients,  but  remaining  near  p  when 
the  gradients  are  small.  Numerical  experiments  with  this  form  of  regulariza¬ 
tion  have  been  very  promising. 

Turbulence  Modeling  for  Optical  Propagation.  This  project  concerns  the 
statistical  characterization  of  atmospheric  turbulence,  from  the  perspective 
of  optical  propagation  through  the  atmosphere.  We  are  currently  collab¬ 
orating  with  PL/LIA  scientists  on  several  topics,  including  the  estimation 
of  mean  and  variance  parameters  for  temperature  fluctuations  in  the  atmo¬ 
sphere,  determination  of  probability  densities  for  these  parameters,  develop¬ 
ing  confidence  bounds  on  parameter  estimates,  and  finding  efficient  simula¬ 
tion  methods. 

The  point  process  models  we  have  developed  for  groundwater  applica¬ 
tions  appear  to  be  applicable  in  atmospheric  turbulence,  in  that  observed 
data  often  exhibit  behavior  that  appears  to  be  nonstationary  in  the  mean. 
Our  MBPP  models  ppear  to  fit  the  data  quite  well,  giving  us  a  reasonable 
decomposition  of  the  time  series  into  regions  of  stationarity. 

Of  crucial  importance  in  the  analysis  of  laser  propagation  is  the  accu¬ 
rate  simulation  of  atmospheric  variability.  In  order  to  perform  these  simu¬ 
lations,  we  need  not  only  estimates  of  parameters  in  the  statistical  model 
of  turbulence,  but  also  an  accurate  sampling  distribution.  We  are  applying 
nonparametic  density  estimation  methods  developed  for  groundwater  studies 
to  the  problem  of  the  sampling  distribution  for  the  phase  screens  (index  of 
refraction  characterizations)  simulated  from  turbulent  atmosphere  statistics. 

Personnel  supported.  Assoicated  with  this  research  effort  are  several 
coworkers.  Dr.  Christian  Wypasek,  a  postdoctoral  researcher,  has  collab¬ 
orated  in  the  statistical  modeling  of  hydraulic  conductivity.  Two  graduate 
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students,  Jeff  Butera  (supported  by  a  DOE  Computational  Sciences  Fel¬ 
lowship)  and  Mike  Jeffris  (both  supported  by  Dept,  of  Education  GAANN 
Fellowships)  have  participated  in  various  aspects  of  parameter  estimation 
in  groundwater  flow  problems.  Butera  has  recently  accepted  a  tenure  track 
appointment  at  High  Point  University,  while  Wypasek  has  been  employed 
since  1  July  97  at  GE  Financial  Services.  Mike  Horton,  a  graduate  student 
in  mathematics  supported  on  this  grant,  is  currently  investigating  on  certain 
statistical  aspects  of  the  ABL  work.  Anand  and  Aravind  Ramachandran, 
two  local  high  school  students  studying  mathematics  at  NCSU,  are  also  pro¬ 
gramming  some  smoothing  algorithms  for  use  on  the  ABL  dataset. 

Publications.  We  have  the  following  publications  to  report: 

•  “Estimation  of  Time  Dependent  Parameters  in  General  Parabolic  Evo¬ 
lution  Systems,”  by  A.S.  Ackleh  and  Ben  G.  Fitzpatrick,  J.  Math. 
Anal.  Appl.  203,  1996,  pp.  464-480. 

•  “Estimation  of  Discontinuous  Parameters  in  General  Nonautonomous 
Parabolic  Systems,”  by  by  A.S.  Ackleh  and  Ben  G.  Fitzpatrick,  Kyber- 
netica,  32,  number  6,  1996,  pp  543-556. 

•  “An  Adaptive  Change  Detection  Scheme  for  a  Nonlinear  Beam  Model,” 
by  M.  A.  Demetriou  and  B.  G.  Fitzpatrick,  Kybemetica,  33,  number 
1,  1997,  pp  103-120. 

•  “On  Approximation  in  Total  Variation  Penalization  for  Image  Recon¬ 
struction  and  Inverse  Problems,”  by  Ben  G.  Fitzpatrick  and  Steven  L. 
Keeling,  accepted,  J.  Num.  Func.  Anal.  Opt. 

•  “Estimation  of  Groundwater  Flow  Parameters  Using  Least  Squares,” 
by  Kendall  R.  Bailey  and  Ben  G.  Fitzpatrick,  accepted,  Math.  Comp. 
Modeling. 

•  “Modeling  Aggregation  and  Growth  Processes  in  an  Algal  Population: 
Analysis  and  Computations”  by  A.  S.  Ackleh  and  B.  G.  Fitzpatrick,  J. 
Math.  Bio.  35,  1997,  pp.  480-502.  Interactions/Transitions. 

•  “On  Continuous  Dependence  under  Approximation  for  Groundwater 
Flow  Models  with  Distributed  and  Pointwise  Observations,”  by  Ben  G. 
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Fitzpatrick  and  Michael  A.  Jeffris,  Discrete  and  Continuous  Dynamical 
Systems ,  2,  no.  1,  1996,  pp.  141-149. 

•  “Sampling  Distribution  of  Approximate  Errors  for  Least  Squares  Iden¬ 
tification,”  by  G.  Yin,  Ben  G.  Fitzpatrick,  and  K.  Yin,  accepted, 
Stochastic  Anal.  Appl. 

•  “Dispersion  Modeling  and  Simulation  in  Subsurface  Contaminant  Trans¬ 
port,”  by  Jeffrey  V.  Butera,  Ben  G.  Fitzpatrick,  and  Christian  J. 
Wypasek,  submitted  to  Math.  Models  and  Methods  in  Appl.  Sci. 

a.  Meetings  and  seminars.  The  Principal  Investigator  has  given  two  talks 
at  Phillips  Lab  worshops  on  Atmospheric  Turbulence  and  Optical  Propaga¬ 
tion  (Dec  96  and  July  97)  and  a  talk  at  the  1997  AFOSR/ARO  Environmen¬ 
tal  Quality  contractors’  meeting.  He  attended  the  Electricity  and  Magnetism 
AFOSR  workshop  at  Brooks  AFB  (Jan  97).  He  has  given  colloquium  talks 
at  Worcester  Polytechnic  Institute  (Mar  97)  and  at  the  University  of  South¬ 
ern  Louisiana  (Apr  97),  and  an  invited  presentation  at  the  1997  Barrett 
Lectures,  University  of  Tennessee  (Mar  97).  He  gave  a  seminar  talk  at  D. 
H.  Wagner,  Associates,  on  the  inverse  interferometry  effort  (July  97).  Dr. 
Christian  Wypasek  and  Mr.  Jeff  Butera  gave  talks  at  the  SIAM  Conference 
on  Geosciences  (June  1997). 

b.  Consultative  functions.  The  Principal  Investigator  visited  Kirtland 
AFB  several  times  during  this  report  period  (Oct,  Dec  96;  April,  May,  June, 
July  97)  to  discuss  various  aspects  of  the  data  analysis  and  modeling  efforts 
for  the  ABL  work.  In  January  of  1997,  he  visited  Brooks  AFB  to  continue 
collaborations  with  scientists  on  E&  M  control  problems.  Several  trips  were 
made  to  Arnold  Engineering  Development  Center  (Nov  96;  May,  July,  Aug 
97)  to  continue  collaborations  on  the  inverse  interferometry  and  design  opti¬ 
mization  efforts. 
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